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Abstract. The thermal Casimir effect in ideal metal rectangular boxes is considered using the method of 
zeta functional regularization. The renormalization procedure is suggested which provides the finite ex- 
pression for the Casimir free energy in any restricted quantization volume. This expression satisfies the 
classical limit at high temperature and leads to zero thermal Casimir force for systems with infinite char- 
acteristic dimensions. In the case of two parallel ideal metal planes the results, as derived previously using 
thermal quantum field theory in Matsubara formulation and other methods, are reproduced starting from 
the obtained expression. It is shown that for rectangular boxes the temperature-dependent contribution 
to the electromagnetic Casimir force can be both positive and negative depending on side lengths. The 
numerical computations of the scalar and electromagnetic Casimir free energy and force are performed for 
cubes. 



1 Introduction 

During the last few years there was an increasing inter- 
est to the Casimir effect pQ, the physical phenomenon 
which is caused by the modification of the spectrum of 
zero-point and thermal oscillations in restricted quantiza- 
tion volumes and in spaces with nontrivial topology (2] 
[3]. This effect has found prospective fundamental appli- 
cations, e.g., for constraining predictions of modern uni- 
fication theories beyond the standard model [4,5,6,7,8,91 
or in the braneworld cosmological scenaria [TUlfTTlfT?] . Al- 
though much of recent research on the Casimir effect was 
concentrated in the experimental issues and material prop- 
erties bi II 1") . some of the principal field-theoretical as- 
pects remain unsettled until the present time. 

One of the most intriguing results is that the Casimir 
energy and force may change sign depending on geometry 
of the configuration and the type of boundary conditions. 
A dramatic example of this situation which has given rise 
to many discussions in the literature for several decades 
is a rectangular box with side lengths a, b and c. Lukosz 
|16j noticed that the electromagnetic Casimir energy in- 
side an ideal metal box may change sign depending on side 
sizes a, b and c. The detailed investigation of the Casimir 
energy for fields of different spins, when it is positive or 
negative, inside a rectangular box as a function of box 
dimensions was performed by Mamayev and Trunov [171 
118) . In particular, the analytical results for two- and three- 
dimensional boxes at zero temperature were obtained by 
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the repeated application of the Abel-Plana formula [2"]ll"3] . 
Ambj0rn and Wolfram [19] used the Epstein zeta function 
to calculate the Casimir energy for a scalar and electro- 
magnetic field in hypcrcuboidal regions in n-dimcnsional 
space-time. The electromagnetic Casimir densities were 
investigated for a wedge [20] and for a wedge with a coax- 
ial cylindrical shell [3T]. The problem of isolation of the 
divergent terms in the vacuum energy and their interpre- 
tation received the most attention. In recent years this 
problem was reformulated as a rectangular box or arbi- 
trary shaped cavity divided into two sections by an ideal 
metal movable partition (piston) 22, 23, 24, 25, 26, 27J. It 
was shown that the Casimir force acting on the piston 
with Dirichlet boundary conditions attracts it to the near- 
est wall. Based on this, some doubts about the previously 
obtained results demonstrating Casimir repulsion in cubes 
have been raised. Both sets of results obtained for stan- 
dard boxes and for boxes with a piston are, however, in 
mutual agreement. The attraction (or repulsion for a pis- 
ton with Neumann boundary conditions [27j ) of a piston 
to the nearest face of the box with any sides a, 6, c does 
not negate the Casimir repulsion for boxes without a pis- 
ton that have some appropriate ratio of a, b and c. The 
point is that the cases with an empty space outside the box 
and with another section of the larger box ouside the pis- 
ton are quite different. In the first case the vacuum energy 
outside the box does not depend on a, b and c and there 
is no force acting on the box from the outside. Whereas in 
the second case there is an extra section of the larger box 
outside the piston which gives rise to the additional force 
acting on it. 

Another intriguing problem to which this paper is de- 
voted is the problem of the thermal Casimir effect in ideal 
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metal rectangular boxes. The hrst calculations on this sub- 
ject |19| resulted in a divergent free energy after removing 
the regularization. More recent results appear to be either 
infinite [28] or ambiguous [29]. Reference [30] reconsid- 
ered the derivation of the Casimir free energy for a mass- 
less scalar and electromagnetic field using zeta functional 
regularization. However, as shown below, the used formal- 
ism does not include all necessary subtractions. This led 
to unjustified conclusions on the behavior of the Casimir 
free energy as a function of temperature. Specifically, it 
is claimed that in rectangular cavities in the high tem- 
perature regime the leading term of the free energy is of 
order (fcgJ 1 ) 4 , where fcs is the Boltzmann constant [30] . 
The terms of orders (AjbT) 3 and (fc^T) 2 are also obtained 
[30 . The coefficients near such terms unavoidably depend 
of a, b and c and, thus, they contribute to the Casimir 
force. 

In the present paper we consider the thermal Casimir 
effect in rectangular boxes starting from the general ex- 
pression for the Casimir free energy in the framework of 
Matsubara quantum field theory at nonzero temperature. 
The thermal correction in this expression is renormalized 
in the same way as the zero-temperature contribution. We 
critically reconsider the known results from a unified point 
of view, separating clearly the empty space contribution 
in the volume of a body and other geometrical terms not 
only in the Casimir energy at zero temperature, but also 
in the temperature correction to it. Using the method of 
zeta functional regularization [31,32,33], the finite expres- 
sion for the Casimir free energy associated with a volume 
V is obtained. This expression is first applied to the fa- 
miliar case of two plane parallel planes and the standard 
results [213,13,34] are reproduced. The Casimir force, as 
defined in this paper, goes to zero when the character- 
istic dimensions of the volume V go to infinity. This is 
achieved by a subtraction of the free energy of the black 
body radiation within the volume V and of two other geo- 
metrical contributions of quantum origin. Then we apply 
obtained general expression for the Casimir free energy to 
the cases of the massless scalar field with Dirichlet bound- 
ary conditions and the electromagnetic field in rectangular 
boxes. We show that the temperature-dependent contribu- 
tion to the Casimir force at zero temperature can be both 
positive and negative (i.e., it leads to repulsion or attrac- 
tion) depending on the box side lengths. This property 
is preserved for boxes with zero electromagnetic Casimir 
force at zero temperature. Numerical computations of the 
scalar and electromagnetic Casimir force as functions of 
side length and temperature for cubes are performed. 

The paper is organized as follows. In Sect. 2 the general 
expression for the Casimir free energy is derived using the 
method of zeta functional regularization. Section 3 con- 
tains the application of this expression to the case of two 
plane parallel planes. In Sects. 4 and 5 the massless scalar 
and the electromagnetic fields, respectively, are considered 
in 3-dimcnsional rectangular boxes. Section 6 contains our 
conclusions and discussion. Throughout the paper the sys- 
tem of units is used where c = K = 1 . 



2 The Casimir energy in Matsubara 
formulation and zeta functional regularization 

In the Matsubara formalism one uses an Euclidean field 
theory, considered as a continuation of the field theory in 
Minkowski space-time by a rotation of the time t — > — it. 
The Euclidean time r is confined to the interval r £ [0, ft], 
where ft — 1/(&bT) is the equivalent dimension corre- 
sponding to the inverse temperature. The bosonic fields 
must be periodic, ip(r + ft, r) — ip(r,r). In the limit of 
zero temperature one re-obtains the theory on the whole 
time axis. 

In the Matsubara formalism, the partition function Z 
has the following representation in terms of a functional 
integral, 

Z = C J D(pe~ SE ^\ (1) 

where Se[<p] is the Euclidean action. It can be obtained 
from the corresponding one in Minkowski space-time by 
the replacement S with iS E - For example, for a scalar field 
with mass m it is 

S E [4>] = lf dT J dr ^Ke <p, (2) 

where 

K E ={-a E + m 2 ). (3) 

The Euclidean wave operator U E = d 2 /dr 2 + A is the 
continuation of the D'Alembert operator. In the functional 
integral |T]), the field to be integrated over must fulfill the 
corresponding periodicity conditions. 

In general, in the Matsubara formalism, the construc- 
tion of the theory, to a large extent, goes in parallel to 
the zero temperature case. This is true for instance in the 
calculation of the functional integral. Since we consider 
free field theories, the functional integral is Gaussian and 
can be calculated directly. Using the infinite dimensional 
analogue to the standard finite dimensional relations, we 
obtain for the partition function from |T]) 

Z = C {dct KeV 1 ' 2 , (4) 

where C is an irrelevant constant which will be dropped 
below. Further, for the free energy we get 

F=~biZ = ±TrlnK E . (5) 

The trace in this expression is taken over the same space 
of fields to be integrated over in the functional integral 
©■ . 

Since we assume thermal equilibrium, it is always pos- 
sible to separate the Euclidean time variable from the spa- 
tial variables. Assuming for the spatial part an eigenfunc- 
tion expansion 

- V 2 <£ 7 (r) = Aj$j(r), Aj = u 2 j-m 2 , (6) 

where J is a collective index, we obtain a basis in the space 
of fields <p, 

$l,j(T,r) = —==$j{r), (7) 

V Z7T 
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with the so-called Matsubara frequencies 
& = 2wk B Tl, I = 0, ±1, ±2, 



(8) 



These functions are periodic in t and these are eigenfunc- 
tions of Ke, 



(9) 



As a consequence, the trace in the free energy becomes a 
sum over the logarithm of the eigenvalues, 



(10) 



=— oo J 



Note that if the field is massless and all quantum numbers 
in the collective index J are discrete, it is assumed in (JTUJ) 
that there are no physical states with Aj = 0. 

In (fTO)) . is the free energy of an ensemble of states 
containing particles at the temperature T. In the special 
case of T — > 0, the time interval stretches over the whole 
axis and the sum over the Matsubara frequencies becomes 
the integral over the frequency £, 



k B T 



OO 

E 



mo 



2tt 



f(0 



(11) 



[here / (£;) is a function which must allow for an analytic 
continuation from discrete values to continuous ones]. In 
this way, J-q defined in (JTUJ) turns into the vacuum energy 

Eq. 

The free energy, as given by (fT0|) , still contains ultravi- 
olet divergences and one has to introduce a regularization. 
In zeta functional regularization, the free energy becomes 

1 r) 00 

^° (s)= 2^ 2s]cbT £ Efe 2 +^+™ 2 r- 



1 — — OQ J 



The regularization is removed for s->0 and fi is an arbi- 
trary parameter with the dimension of mass. The separa- 
tion of the ultraviolet divergences can be done quite easily 
because these are the same as at zero temperature. There 
are two ways to proceed. In the first method one has to 
apply the Abel-Plana formula 2, 13 to the frequency sum 
in (fT2)l . This way has the advantage that it can also be ap- 
plied to the case when there is an additional dependence 
on I, for example through the dielectric permittivity enter- 
ing in the form of £ 2 — > e(i£)£ 2 - Another, to some extent 
easier, way is through the application of the Poisson sum- 
mation formula [35] . According to this formula, if c(a) is 
the Fourier transform of a function b(x), 



c(a) 



then it follows that 



1 

2^ 



b(x)e 



c dx, 



£ 6(0 = 2tt £ c(2nl). 



(13) 



(14) 



By putting 
b(x)-- 



c(a) 



1 



2^/ttz 

we obtain from (fT4")l the following equality: 



E< 

1 — — OQ 



,-zr 



7 E 



(15) 



(16) 



n— — oo 



where Re z > is assumed. In order to use this equality 
we represent (|12p as a parametric integral, 



j2 ^ e - t («' +Aj+m2 ) 



1 — —OQ J 



™ --was" I TW) 

(17) 

and apply {TSJ with z = (2nk B T) 2 t = {2n/(5) 2 t. The re- 
sult is 



Ms) 



ld_ 
28^' 



,2s 



E 



n— — oo 



dt 



t r(s)V47rt 



E< 



x 



jf t(Aj+m 2 ) 



(18) 



The n-dependent factor in the exponential provides 
convergence for the ^-integration at t — > for all terms in 
the sum over n except for n = 0. The latter is just the zero 
temperature contribution. This can be seen by applying 
(fTTj) to the frequency sum in (fTT)) . 



k B T £ 



2tt 



1 



/Ant 



(19) 



As a consequence, we can split the free energy into the zero 
temperature part and the temperature dependent addition 
AtTq as 

F {s) = £J ,eff(s) + A T T {s). (20) 
Here, the vacuum energy at zero temperature is 



The nonrenormalized thermal correction is given by 



(21) 



0_ 
d~s 



dt 



t a 



t r(«)V47rf 

^ e -^ 2 -*( Aj+m2 ). (22) 



/ — — oo 



Note that (|20[) has a transparent physical interpretation 
only for the temperature independent boundary condi- 
tions considered here. 

The ultraviolet divergences are contained in Exj, e ff(s) 
and can be dealt with in the standard way [2lll3j. This 
results in the replacement of Eo >e g(s) with i?J en in (|2"0|) . 
where E^ en is the finite renormalized Casimir energy at 
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zero temperature associated with a volume V, which goes 
to zero when V — ► oo. In AtJ 7 o(s) in (|22]) . the integration 
over t is convergent and we can remove the regularization, 
i.e., we can put s = using the equality 



m 

6 ds r(s) 



lim — 



/(0), 



(23) 



where /(s) is any regular function at s — 0. Following this, 
the integration over t and the summation over n can be 
carried out explicitly 



A T T Q = k B T^\n(l- e 



(24) 



In this formula, we see the sum of the T-dependent con- 
tributions to the free energies of the individual degrees of 
freedom, or modes Aj, of the considered system. Taking 
into account Aj + m 2 = uij, the total free energy of all 
the oscillator modes appears to be 



•^0 — 



fc B rVln(l -e 



j 



where the zero temperature contribution is already re- 
placed with Eg en . For instance, if we take the volume V 
to be that of empty space, the modes are plane waves, the 
index J becomes the wave vector k and the sum over J 
turns into the corresponding momentum integration with 
respect to dfc/(27r) 3 . As a result, from (j2~4")l we obtain the 
free energy density of the blackbody radiation 



/bb(r) = *BT/ * ln(l-e-^l 



7T 2 (fc B T) 4 



90 

(26) 

We note that (|26f holds for a scalar field. For the electro- 
magnetic field one would have to multiply by a factor of 
2. In fact, for the empty space /bb defined in (|26|) is the 
complete free energy. This is because in this case the zero 
temperature part was the vacuum energy of empty space 
which was disregarded in the replacement of Eo. g with 
£^ cn . Using the thermodynamic connection between the 
energy at a temperature T and the free energy 



U(T) = -T' — 



2 d T{T) 



T 



(27) 



it is evident that for T{T) = f^{T) = 2/ bb (T) the 
respective energy density is in agreement with Planck's 
blackbody radiation density 



ir 2 (k B T) 4 
15 



(28) 



If we consider the free energy in a restricted volume 
V, then in accordance with (|25[) . we have to keep the zero 
temperature part E^, where the empty space contribu- 
tion has been dropped already. For the temperature de- 
pendent part, we have to note that we are interested in 
the change in energy which comes from the volume V . 
Therefore we need to subtract from the temperature de- 
pendent part AtTq of the free energy, the corresponding 



amount related to empty space, i.e., the blackbody radia- 
tion density /bb multiplied by the volume V. As a result, 
we come to the following expression for the renormalizcd 
free energy associated with a finite volume V, 



T 



+ A T F - Vh 



bb. 



(29) 



which can be considered as the relevant quantity for the 
Casimir effect at finite temperature. In the next section 
we show that this is the case for the configuration of two 
parallel ideal metal planes. 

However, in general case (|29|) is not satisfactory and 
requires at least two additional subtractions in order to 
get the physical Casimir free energy. The point is that the 
asymptotic expression of the thermal correction, AtJ-o, 
in the limit of high temperatures (large separations) in 
addition to the leading term, 



V 



7T 2 (k B T) 4 



90(ftc) 



(25) contains two other similar terms 



{ k B Tf 
{Tic) 2 



OL 2 - 



(k B T) 2 



he 



(30) 



(31) 



where a\ and a 2 have dimensions of cm 2 and cm, re- 
spectively. They are expressed through the heat kernel 
coefficients a\/ 2 and a\, respectively |13j . This was shown 
quite generally by Dowker and Kennedy [36] . We use usual 
units in (|3"D)) and (|3"Tj) to underline the quantum character 
of all these terms. The coefficients ct\ and generally 
speaking, depend on geometrical parameters of the con- 
figuration (for two parallel planes a\ = a 2 = 0, but for a 
rectangular box a.\ and a 2 depend on the side sizes; see 
Sects. 4 and 5). Because of this, the terms of the form 
(f3"Tj) in the free energy would lead to forces of quantum 
nature which do not vanish with increasing characteris- 
tic sizes of the body (the next expansion term in the 
high-temperature limit of the free energy has the form 
of a^k B T with a dimensionless coefficient a^; it is of clas- 
sical origin and does not contribute to the Casimir force). 
Below we demonstrate that for rectangular boxes the ge- 
ometrical structure of the coefficients ai and a 2 is just 
the same as in the two infinite terms subtracted from the 
zero-temperature Casimir energy to make it finite in the 
cutoff regularization |2,13J. Thus, one can perform addi- 
tional, finite, renormalization of the free energy resulting 



^phys = ^rcn + ^ (3^ 

A T T = AtJ-q - Vf hh - ax{k B Tf - a 2 (k B T) 2 . 



Here, the quantity A^T is the physical thermal correction 
to the Casimir energy. The respective Casimir force obvi- 
ously vanishes when the characteristic sizes of the volume 
V along all three coordinate axes go to infinity. 

Below we use (1321) to investigate the thermal Casimir 
effect in configurations of two parallel ideal metal planes 
(Sect. 3) and rectangular boxes (Sects. 4 and 5). All other 
thermodynamic quantities like force and entropy can be 
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also derived from this formula. As for configurations con- 
taining translational invariant directions like parallel planes 
or a cylinder, one should bear in mind that all quantities 
in (|32p must be divided by corresponding parameters like 
the area of plates. 



3 The electromagnetic Casimir free energy 
between two parallel planes 

The Casimir free energy in the configuration of two par- 
allel ideal metal planes can be easily obtained from the 
general equation (|32p . It is common knowledge that for 
two parallel planes the high-temperature asymptotic ex- 
pansion of AtJ~o does not contain terms presented in (|3ip . 
Then (|32|) and (|24| with oi\ — a 2 — represent the elec- 
tromagnetic Casimir free energy per unit area in the form 
(prime adds a multiple 1/2 to the term with n = 0) 



F^ ys (a,T) = - 



720a 3 



(33) 



1 rri pOC 00 . 



-Pu k 



n 2 (k B T) 4 
45 



where 



By introducing the new variable 

z = j3uj k . n = 



1-uT » ^ 



Trn\' 



(34) 



(35) 



and changing the order of summation and integration, we 
rearrange ([3"3"]) as 



F^ ys (a,T) 



+ 



(k B T) 



3 00 , 

■E 

n=0 



720a 3 

z dz ln(l 



(36) 



') + 



720aH 4 '■ 



where t = T e s/T, and the effective temperature is defined 
in the usual units as k B T e s = he /(2a). The integral en- 
tering (|36p can be evaluated using the series expansion 



: dz ln(l 



OO 

1=1 



z dz e 



-iz 



8na 



-2-rrnlt 



(37) 



i=i 



Substituting this into f|36|) and performing the summation 
in n, we obtain 



^ ys (a,T) 



720a 3 



/=1 



3th(7T/t) 



(38) 



t 2 Z 2 sinli (TrfZ) 



This is a well known result [2,3, 13,3407], and the last 
term on the right-hand side of (|38p originates from the 
subtracted contribution of blackbody radiation. At low 
temperature T < T cS (|3gJ> leads to [21 [3TrT3TI54l [37] 



^S hys (a,T) 



720a 3 



45C(3) 



T 



T 

T c s 



(39) 

where £(z) is the Riemann zeta function, and exponen- 
tially small terms are omitted. We emphasize that the 
quantity l/(aT e ff) 3 does not depend on a. Because of this, 
only the last term on the right-hand side of (|3"9")) . origi- 
nating from the subtraction of the blackbody radiation in 
(f32|) . contributes to the measurable quantity, the Casimir 
pressure following from (|39|) 



P e i(a,T) 



df£ ys (a,T) 



1 



T 

T c s 



da 240a 4 

'(40) 

Thus, if this subtraction would not be done in (|32p and 
(|33p . only the exponentially small thermal corrections of 
order (T/T c s) exp(— 2ttT c s/T) were added to unity on the 
right-hand side of (l40|) . 

In a similar way, at high temperature T 3> T e s 
results in [2l3lfT3ll34ll37] 



k R T 



^ e r s (a,T) = -^C(3), Pei(a,T) = 



k B T 
47ra 3 



C(3). 



( 41 ) 

This is the commonly known classical limit where the 
result does not depend on the Planck constant [38ll39j . 
It is obtained because at high temperature the free en- 
ergy of the blackbody radiation is canceled by the other 
terms contained in ([38]) . If the subtraction of the black- 
body radiation in ([32]) were not done, the term of order 
(k B T) 4 /(ch) 3 (in usual units) would be dominating at 
high temperature (or large separations) in contradiction 
with the classical limit. 



4 Scalar Casimir effect in rectangular boxes 

After the above demonstration that the general equation 
([32l) leads to the familiar result in the case of two parallel 
planes, we apply it to calculate the thermal Casimir effect 
in rectangular boxes with sides a, b and c. We start with 
the case of a massless scalar field with Dirichlet boundary 
conditions. Then (13"21), (1241) can be written in the form 



T phys (a,b, c, T) = E 1 cn (a,b, c) 

OO 

+ k B T - e- 0Unl ") 

n,/,p— 1 

^(k B Tf 



!)() 



abc- ai (k B TY -a 2 (k B TY, (42) 



where 



^lip^n 2 [— + — + n, I, p= 1, 2, 3, ... (43) 
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and the Casimir energy at zero temperature is given by 
[40,41, 



E™(a,b,c) 



7T 2 foc 

1440a 3 



C(3)(6 + c) 
327ra 2 



7T 

96a" 



7T 

2a 



a \ a a 



rw +c S 

Here the following notations are introduced 

1 oo 

oo 

R(z 1 ,z 2 ) = -g- E (1 - i5 ;o 5po) 

— OO 

=o / . \ 3 / 2 , 



(44) 



(45) 



and K n (z) are Bessel functions of imaginary argument. 

Equation (|42|) presents the finite renormalized value 
of the Casimir free energy for a scalar field with Dirichlet 
boundary conditions in a box with sides a, b and c valid at 
any temperature. However, the values of a.\ and a 2 remain 
unknown. To determine them we should find the asymp- 
totic expression for the nonrenormalized thermal correc- 
tion At To at high temperatures (large separations). To 
do this, we identically rearrange AtTq to the form 

A T F (a, b, c, T) = k B TX(p a , ft, ft), (46) 

OO 

X(ft,ft,ft) = £ In(l-e-vWW), 

r 

where 



n,J,p— 1 



fl - ^ fl - ^ Ft - ^ 
Pa — 1 P6 — -7-) Pc — ■ 

a b c 



Note that the quantity X does not depend on a, 6, c and 
T separately, but only through the products aT, bT and 
cT. Below we find the asymptotic expression of X under 
the conditions ft, /3ft, ft -C 1. This can be done by the 
repeated application of the Abel-Plana formula [2"1IT3] 



E/w = -o/(°)+ / / 

n=l 2 - / ° - 70 



First we put 

oo 

/(n) = £ In (l - e~V?> 2 +W 2 +^ . 
l,p=i 

Then the application of (|48|) leads to 

+ J dt £ In (l- e-V^ 2 +^ 2 +^ 



(48) 
(49) 



i,p=i 

+0(lnft,lnft,lnft). 



(50) 



Here, it is taken into account that the last term on the 
right-hand side of (|4"B")) with / defined in (|4"9"| is of order 
In ft, In ft and In ft. 

Applying the Abel-Plana formula to each of the sums 
in (fSTJ)) , we get 



oo 

X((3 a ,p b ,p c ) = ~^ln(l-e-^) 

i />00 /» CO °° 

I dy I dv Vlnfl 



(51) 



ft 

+ 0(lnft,lnft,lnft). 
Bearing in mind that 



£ In (1 - e"^) = + 0(ln/3 c ) (52) 

P =i c 

and applying the Abel-Plana formula to the remaining two 
sums, the following result is obtained: 



X(ft,ft,ft) 



24ft 



(53) 



-(— — 

4 \P a + fty jo 

1 / 1 1 

2 Ua/^c + ftft + ft Pb 



(i-e-y) 



dy In (1 — o 
1 



dy 



dv In 1 



(> 



-| />00 /"OO /<OG 

<«» + Kml dy l dv l d " bl ( 1 - e ~' / ^) 



+ 0(lnft,lnft,lnft). 

Calculating all integrals and using (|4*6"|) and notations (|4"T)) 
we arrive at the asymptotic expression for the nonrenor- 
malized thermal correction at high temperatures (large 
separations) 



A T To{a,b, c, T) 
C(3) 



--(fc B T) 2 (a + 6 + c) (54) 



+ 



(ac + bc + ab)(k B TY 



4tt /x ' 90 

0(fcijT In ft, fc B T In ft, £; B T In ft) 



(ksT) abc 



Thus, we have demonstrated that at high temperatures 
(large separations) the nonrenormalized thermal correc- 
tion really contains the terms of form (|30[) and (I?!]) . From 
the comparison of (f54|) with (f3Tj) , it follows 



C(3) 

47T 



(flc + foe + afo), a 2 



24 



(a + fo + c), (55) 



i.e., ol\ is proportional to the total area of box surface, 
and ct2 to the sum of the sides. These geometrical ob- 
jects (together with a box volume) are usually renormal- 
ized when replacing E , c g with E r en [2 | ll3 1 ll7 ( ll8j . Thus the 
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subtraction of the last three terms on the right-hand side 
of (|42|) can be interpreted as an additional finite renor- 
malization giving a physically meaningful temperature- 
dependent contribution to the Casimir energy. 

As mentioned in Sect. 2, ol\ and a 2 can be expressed 
in terms of the heat kernel coefficients. Thus, keeping in 
mind that the heat kernel coefficient ai/ 2 = — ^/ttS/2 [42] , 
where S is the area of box surface, and comparing (I54|) 
with the general asymptotic expression for the free energy 
|36j , one arrives at 



C(3) 

1/2 



«3) 
8tt 



S 



(56) 



in agreement with (|55|) . Similar comparison of the general 
expression [36] with (|54| leads to a 2 = — di/24. The heat 
kernel coefficient a\ can be calculated from the known 
expression for the heat kernel coefficient of an angle 9 43] 



d(0) 



e 2 



(57) 



In the case of rectangular box = ir/2 and c\ = 7r/4. In 
the 3-dimensional case this should be multiplied by the 
length of all sides leading to 



ai = 4ci(a + b + c) = ir(a + b + c). 



(58) 



From this, the expression for a 2 in (|55[) is reobtained. 

The Casimir force acting between the opposite faces of 
the box is obtained from (|4*2")l and (|55|) 

F x (a, b, c,T) = - ^ PhyS( °' b ' C ' r) =F x (a, b, c) (59) 



da 



- y 



7T 2 (k B TY 



90 



be 



+ M {kB Tf{b + c)-^{k B T)\ 

It is well known [T3TI4T] that the scalar Casimir energy 
i?Q Cn (a, b, c) is negative and the respective force F x (a, b, c) = 
—dE(a,b,c)/da is attractive for any ratio of a, b and c. 
Because of this, we restrict ourselves to the consideration 
of a cube a = b = c. 

In the limit of low temperature T <C T e g the leading 
terms in (l42l) are 



F phys {a, T) = E™ n {a) - k B T exp 
, ^{k B Tf 3 3C(3) 



90 



" -—(k B Tfai + -(k B Tra, 



7T\/3 v 

ak B T 1 
n 

2A' 



(60) 



where the scalar Casimir energy for a cube can be calcu- 
lated numerically using (|44p and (|45|) with the result 



0.0102 



(61) 



At arbitrary temperature it is convenient to represent 
both the free energy and force in terms of the dimension- 
less variable t 

F phys (a,T)=E r Q cn (a) 



+ _L m(l-e-W« 2 +< 2 +P 2 ) 



n,l,p—l 



7T 2 3C(3) 1 7T 1 

1440a< 4 ~ 32tt ^ + 32 ^ : 



F x {a,T) = F x {a) 



(62) 



+ 



,2 



-^ =1 ^n 2 + l 2 +p 2 e 2^V« 2 +i 2 +P 2 _ i 



C(3) 1 



7T 1 



1440a 2 t 4 16tt a 2 t 3 96 a 2 t 2 ' 

where the force at zero temperature in accordance with 
(fBTj) is given by 



F x (a) 



0.0102 
3a 2 ' 



(63) 



In Fig. 1(a) we plot the scalar Casimir free energy in 
a cube from (|6"2")l as a function of the side length a at 
T = 300 K (solid line). As is seen from this figure, the free 
energy increases monotonously with the increase of a. At 
large separations (not shown in the figure) it approaches to 
a constant. In the same figure the Casimir energy Eq cu = 
EQ Cn (a) at zero temperature is shown by the dashed line. 
In Fig. 1(b) the scalar Casimir free energy is plotted as a 
function of temperature in the cube with a = 2 /xm. It is 
seen that at large temperatures :F phys is proportional to 
the temperature in accordance with the classical limit. 

The magnitude of the scalar Casimir force from (f62|) 
(in a logarithmic scale) as a function of a at T — 300 K is 
shown in Fig. 2(a) by the solid line. The force is attractive 
for cubes with any side length and its magnitude goes to 
zero with the increase of a. In the same figure the dashed 
line shows the magnitude of the Casimir force (in a loga- 
rithmic scale) acting on the opposite faces of the cube at 
T = 0. Fig 2(b) shows the Casimir force as a function of 
temperature for the cube with a = 2 /an. It is seen that 
for both negative and positive values of the free energy 
the respective Casimir force is attractive. 



5 Electromagnetic Casimir effect in 
rectangular boxes 

Now we consider the electromagnetic thermal Casimir ef- 
fect in a rectangular box with the side lengths a, b and c. 
For the electromagnetic field the renormalizcd free energy 
(|32|). ([24)1 is specified as 



Ff* s {a,b,c,T) = El%a,b,c) 



(64) 



k n T 



E M 1 



-Pu„ 



l,p=l 



„Z=1 



J2 H 1 - e " /3 "" p ) + 2 E - e~^) 

n,p—l n,l,p—l 
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7T 2 (k B TY 

45 



abc - af{k B Tf - af{k B Tf 



Here, uo n i p is defined in (j4"3")l , ui n i — oj n iOi ancl a i \ a f have 
to be determined. The electromagnetic Casimir energy at 
T = 0, E r cn {a, b, c), is given by 



35S(M,c) 



ir 2 bc 



C(3)C 7T /l 



720a 3 16tt6 2 48 



a a 



(65) 



where the functions G and i? are defined in (|45|) . 

In order to find the coefficients af and af, one should 
determine the asymptotic behavior of AtFq^i in the limit 
of high temperatures (large separations). In the electro- 
magnetic case the nonrenormalized thermal correction can 
be identically rearranged to the form 

A T f Q<A (a, b, c, T) = k B TY{(3 a , (3 b , #.), (66) 
Y(p a ,0 b ,0 D ) = 2X(p a ,0 b ,0c) 

+ £ In (l - e-VWTM^ 
i,p=i 

+ f; In (l - c-v^™ 2 +^ 2 ) 

oo 

+ £ In (f - e -V/3> 2 +^V) 5 

n,p— 1 

where the asymptotic behavior of X((3 a , /3 C ) at small 
/3 a , /3ft and /3 C was already determined in Sect. 4. Taking 
into account (|46]) . (|47l) and (I54"]) , it is given by 



, <(3) 



24 V/3 a A fl 
1 1 1 



4 VAA /?a/3c PbPc 

r 5 f 



90 



+ 0(ln/3 Q ,ln/3 6 ,ln/3 c ). (67) 



The asymptotic behavior of Y((3 a , (3 bl Pc) at small f3 a ,fl b , Pc 
is obtained in perfect analogy with the case of the scalar 
field by the repeated application of the Abel-Plana for- 
mula (|48[) to the remaining three summations in (|66p . The 
result obtained with account of (|6"T)) is 



Y(p a ,0b,0 c ) 



1 



45 



12 \p a (3 b /3 C 



+ 0(ln/? a ,ln/3 h ,ln/? c ). 



(68) 



Substituting this into ([66]) and using notations (|47l) . one 
obtains the asymptotic expression for the thermal correc- 
tion at high temperatures (large separations) 

A T F 0M (a, b, c, T) = ^{k B T) 2 {a + b + c) (69) 



45 



(k B T) 4 abc+0(k B Tlnl3 a ,k B Tln/3 b ,k B Tlnl3 c ). 



It is notable that |69|) does not contain a contribution 
proportional to the surface area of the box [compare with 
(|54p for the scalar field] . In the electromagnetic case such a 
contribution is absent also in the divergent Casimir energy 
of the box Eo >e i at zero temperature [2l H31[T71fT5] . Thus, 
from (l69l and (I3T1) one arrives at 



af = 0, 



12 



(a + b + c). 



(70) 



The Casimir force acting between the opposite faces of 
a box is obtained as the negative derivative of (|64|) with 
respect to a, 



F x ,ei(a, 6, c, T) = i^ei (a, 6, c) + 



, , w„z(e- 

n,t— 1 V 



„, p =i w «p 

n 2 (k B TY 
45 



\ " 



- 1) 



n,l,p— 1 



- =1 w„ ip (e^"'p 



6c- 



(71) 



It has been known that the electromagnetic Casimir 
energy inside a box at T — 0, E^(a,b,c), can be both 
positive and negative and the Casimir force, F XiC i(a, b, c) = 
—dE^ A {a,b,c)/da, can be both attractive and repulsive 
depending on the ratio of the sides a, b and c [21IT3]. Here 
we consider in more detail the thermal electromagnetic 
Casimir effect for a cube a = b = c where the electro- 
magnetic Casimir energy at zero temperature computed 
by (|65|) is positive, 



(72) 



and the force (|7Tj) is repulsive. 

For a cube the electromagnetic Casimir free energy 
(f6"4"| and force (fTTj) are given by 



^(a,T)=E^(a) + ^ t J2 ln(l 



-2-7TtVn 2 +l 2 



n,l=l 



n,l,p— 1 



F x ,ei(a,T) = F S) ei(a) 



(73) 



2vr 

7^2" 



E 



1 



n,l=l 



Vn 2 + Z 2 e 27r^ ^ /7;T + ^2 " - 1 



720a 2 i 4 48a 2 t 2 ' 
where the force at T = is 

Fx,el( a ) ' 



0.09166 
3a 2 ' 



(74) 
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In Fig. 3(a) we plot the electromagnetic Casimir free 
energy in a cube as a function of a at T = 300 K (solid 
line). In the same figure the Casimir energy at T = 
is shown by the dashed line. As is seen in this figure, the 
electromagnetic Casimir free energy decreases with the in- 
crease of separation. Similar to the scalar case, at large 
separations ^"P hys approaches to a constant. In Fig. 3(b) 
the electromagnetic Casimir free energy is shown as a 
function of temperature for a cube with a = 2 /jm. The free 
energy decreases with the increase of T. At high temper- 
atures .Fgj hy8 demonstrates the classical limit. The respec- 
tive thermal electromagnetic Casimir force at T = 300 K, 
as a function of a, is shown in Fig. 4(a) by the solid line. It 
is positive (i.e., repulsive) for cubes with any side lengths. 
Thus, thermal effects for cubes in the electromagnetic case 
increase the strength of the Casimir repulsion. The dashed 
line in Fig. 4(a) shows the electromagnetic Casimir force 
at T = as a function of a. This force is given by (|74| . i.e., 
it is always repulsive. Fig. 4(b) demonstrates the electro- 
magnetic Casimir force in a cube of side length a = 2 fim 
as a function of temperature. It is seen that the force in- 
creases with increasing temperature. 

Note that the obtained results differ from those found 
in [3U], where the terms of order (fc^T) 4 and of lower or- 
ders in the Casimir free energy were obtained in the high- 
temperature regime. Also, the Casimir free energy in [30) 
both for the scalar and electromagnetic field is always a 
decreasing function of temperature on the opposite to our 
result in Fig. 1(b). This is explained by the fact that in 
[3"0] the subtractions of the contributions from the black- 
body radiation and of the terms proportional to the box 
surface area and to the sum of its sides were not made. 

The above equations ([42jl. ([59]) . (|64j) and ((TTJ) can be 
used to compute the scalar and electromagnetic free en- 
ergy and force for boxes with arbitrary side lengths a, b 
and c. Specifically it follows that the temperature-depen- 
dent contribution to the electromagnetic Casimir force 
[which is obtained as —dAxT/da from the physical ther- 
mal correction defined in l|32p] can be both positive and 
negative depending on the side lengths a, b and c. On the 
one hand, as shown above (see Fig. 4), for a cube ax ax a 
the temperature-dependent contribution to the Casimir 
force is positive and computations show that this is pre- 
served for any box a x b x b with a > b. On the other 
hand, for boxes with b = c = 10 /xm and a\ — 2.942 /im 
or d2 = 34.29 ^m the Casimir energy at T = is equal 
to zero [2lfT3], Computations using l[7T|) show that for the 
box a\xb xb the temperature-dependent contribution to 
the force is negative, whereas for the box ai x b x b it is 
positive. 

The thermal correction to the Casimir energy and force 
acting on a piston were investigated [33] for the scalar 
field with Dirichlet or Neumann boundary conditions us- 
ing (25). The electromagnetic Casimir free energy and 
force acting on a piston were found in the case of ideal 
metal rectangular boxes and cavities with general cross 
section 24 . In the limit of low temperatures the thermal 
correction to the Casimir force on a piston was shown to 
be exponentially small. In the limit of medium tempera- 



ture [24] obtains the term of order (fc^T) 4 and of lower 
orders in the electromagnetic Casimir free energy (see (44) 
and (47) in [24]). This results in the contribution to the 
force acting on a piston which increases with the increase 
of the temperature, depends on Ti and c and does not de- 
pend on the position of a piston. The same results are 
obtained from the application of ([3"2"|) . This is because the 
contribution of the blackbody radiation to the energy of 
entire box is equal to — a&c/bb and does not depend on 
the position of a piston. The contributions from the terms 
proportional to the surface areas and the sums of the side 
lengths for the two boxes to the left and to the right of a 
piston are also independent of the piston position. 

6 Conclusions and discussion 

In the above we have reconsidered the old problem of 
thermal Casimir force in rectangular boxes made of ideal 
metal. Our main concern was with the appropriate defini- 
tion of the Casimir free energy inside some closed volume. 
We have stressed that such a definition must include a 
subtraction of the energy of blackbody radiation confined 
in this volume and of two other terms proportional to 
(fcsT) 3 and (ksT) 2 which are contained in the high tem- 
perature (large separation) asymptotic expression for the 
thermal correction to the Casimir energy AtTq in (|2"4")1 . 
If this is done, the thermal Casimir force turns into zero 
when the characteristic sizes of the system go to infinity. 
According to the results of Sect. 3, the suggested expres- 
sion (f3"2"| for the Casimir free energy is applicable to the 
limiting case of two plane parallel ideal metal planes and 
leads to the same result as obtained by many authors in 
the framework of thermal quantum field theory in Mat- 
subara formulation. Both the free energy and the thermal 
Casimir force were calculated for the massless scalar field 
and for the electromagnetic field inside ideal metal rectan- 
gular boxes. For the scalar field with Dirichlet boundary 
conditions the thermal Casimir force is always attractive. 
It was shown that for the electromagnetic field in rect- 
angular boxes the temperature-dependent contribution to 
Casimir force can be both positive and negative depend- 
ing on the side lengths. At high temperatures the Casimir 
free energy is proportional to fc^T, i.e., the classical limit 
is achieved, like it holds for two plane parallel planes. The 
suggested expression for the Casimir free energy leads to 
the cancellation of all terms of order (fc^T) 4 , (fcsT) 3 and 
(ksT) 2 at high temperatures (large separations). This is 
in accordance with what the physical intuition suggests 
because such terms depend on the Planck constant and, 
thus, are of quantum origin. 

The above results are obtained for rectangular boxes 
with the Dirichlet boundary conditions (scalar case) or 
for ideal metal boxes (electromagnetic case). We would 
like to emphasize that the final resolution of the problems 
of Casimir repulsion and of the thermal Casimir effect in 
rectangular boxes requires the inclusion of real material 
properties of the boundary surfaces. In the case of closed 
quantization volumes this is a challenging issue which re- 
mains to be explored. 
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Fig. 1. The scalar Casimir free energy for a cube as a function of (a) a side length a at T = 300 K (solid line; the dashed line 
shows the energy at T = 0) and (b) temperature at a — 2 fim. 
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Fig. 2. The scalar Casimir force between the opposite faces of a cube as a function of (a) a side length a at T = 300 K (solid 
line; the dashed line shows the force at T = 0) and (b) temperature at a = 2 fim. 
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Fig. 3. The electromagnetic Casimir free energy for a cube as a function of (a) a side length a at T = 300 K (solid line; the 
dashed line shows the energy at T = 0) and (b) temperature at a — 2 /im. 
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Fig. 4. The electromagnetic Casimir force between the opposite faces of a cube as a function of (a) a side length a at T = 300 K 
(solid line; the dashed line shows the force at T = 0) and (b) temperature at a — 2 fim. 



